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A Class  of  Totally  Positive  Blending  B-Bases 


Laura  Gori,  Laura  Pezza,  and  Francesca  Pitolli 


Abstract.  Totally  positive  blending  bases  present  good  shape  preserving 
properties  when  they  are  used  in  CAGD.  Among  these  bases  there  exist 
special  bases,  called  B-bases,  which  have  optimal  shape  preserving  prop- 
erties. In  particular,  the  corresponding  control  polygon  is  nearest  to  the 
curve  among  all  the  control  polygons;  thus  many  geometrical  properties 
are  similar  to  the  ones  of  the  curve.  Examples  of  totally  positive  blend- 
ing B-bases  are  the  Bernstein  polynomials  and  the  B-spline  basis.  Our 
purpose  is  to  construct  new  classes  of  such  bases  starting  from  compactly 
supported  totally  positive  scaling  functions. 


§1.  Introduction 

One  of  the  main  goals  in  Computer  Aided  Geometric  Design  (CAGD)  is  to 
predict  or  control  the  shape  of  a curve  by  studying  or  specifying  the  shape  of  the 
control  polygonal  arc  formed  by  certain  points  which  define  the  curve,  typically 
the  coefficients  when  the  curve  is  expressed  in  terms  of  a particular  basis.  This 
is  possible  when  we  choose  as  a basis  a system  of  functions  v = («o>  • • ■ , vn) 
with  suitable  shape  preserving  properties.  This  means  that  the  geometrical 
properties  of  the  curve  in  1R2 


7(z)  = ^ PiVi(x),  x € I C 1R,  (1) 

i=0 

constructed  on  the  control  points  Pi  € 1R2,  i = 0, . . . ,n,  are  implied  by  the 
geometrical  properties  of  the  control  polygon  P0...Pn.  The  shape  preserving 
properties  of  each  representation  (1)  depend  on  the  characteristic  of  the  system 
v. 

The  bases  commonly  used  in  CAGD,  such  as  Bernstein  bases,  B-splines, 
/3-splines,  nonuniform  rational  splines  (NURBS),  are  blending  totally  positive 
systems.  This  means  that  the  collocation  matrix 
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for  any  sequence  xo  < ■ • • < xs,  Xi  6 I,  i = 0, . . . ,s,  is  totally  positive  (i.e., 
all  its  minors  are  non- negative),  and  the  basis  functions  add  to  one,  that  is 


n 

]r\i(z)  = l,  x e 7.  (3) 

i=0 

The  importance  of  blending  totally  positive  systems  is  due  to  the  fact  that 
they  enjoy  two  properties  which  are  usually  demanded  for  curve  control:  the 
convex  hull  (CH)  and  the  variation  diminishing  (VD)  properties  (see,  for  in- 
stance, [5,6]).  As  a consequence,  in  many  ways  the  shape  of  the  curve  7 
mimics  the  shape  of  the  control  polygon  P0  . . . Pn.  However,  blending  totally 
systems  usually  do  not  enjoy  a property  which  is  also  important:  the  end-point 
interpolation  (EPI)  property. 

Bases  which  simultaneously  satisfy  the  VD,  CH  and  EPI  properties  can 
be  obtained  by  considering  blending  B- bases  [5]. 

Following  [5],  a totally  positive  (TP)  system  u of  linearly  independent 
functions  is  said  to  be  a B-basis  if  any  totally  positive  basis  v of  the  space  U 
generated  by  u satisfies  the  condition 

v = uA,  A nonsingular  totally  positive  matrix.  (4) 

In  [4]  it  was  proved  that  if  there  exists  a blending  TP  basis  in  U,  then  there  ex- 
ists a unique  blending  B-basis  for  that  space.  B-bases  have  optimal  properties 
in  the  geometric  context  [5],  that  is,  in  particular,  the  control  polygon  with 
respect  to  the  B-basis  is  nearest  to  the  curve  among  all  the  control  polygons 
with  respect  to  any  other  TP  basis. 

Some  examples  of  B-bases  are  given  in  [4,5];  in  particular,  the  B-spline 
basis  is  the  blending  B-basis  in  the  space  of  the  polynomial  splines  of  degree 
m,  on  a given  interval  with  a prescribed  sequence  of  knots. 

At  this  point,  it  is  worthwhile  to  remark  that  in  the  case  of  cardinal 
splines  (knots  at  the  integers),  this  basis  is  related  to  the  cardinal  B-spline 
Nm,  defined  by  Nm  = Nm~l  * N°,  where  N°  is  the  characteristic  function  of 
[0, 1)  and  * denotes  the  convolution  product  (see,  for  instance  [8]). 

On  the  other  hand,  Nm  is  a scaling  function,  that  is  the  solution  of  the 
functional  equation 


- 771-i-l  / I 1 \ 

Nm(x)  = - )^m(2x-f),  *€R.  (5) 

i= 0 ' ' 


In  this  paper,  we  analyse  the  more  general  problem  of  the  construction 
of  blending  B-bases  considering,  instead  of  Nm,  a scaling  function  satisfying 
a functional  equation  more  general  than  (5): 

p(x)  = XX2*-*-), 

i£2Zi 


X E R, 


(6) 
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a2i+l  = 52  ~ (7) 

ieE  ieZ 

It  is  known  that  a solution  ip  of  (6)  exists  if  the  mask  a satisfies  further 
conditions,  in  addition  to  (7).  In  particular,  if: 
i)  a is  compactly  supported  on  [0,  m + 1]  (with  aoam+i  / 0), 
ii)  the  symbol 

771+1 

P(z ) = 52  aizi  (8) 

t=0 

has  roots  with  negative  real  part  (Hurwitz  polynomial), 
then  there  exists  [8]  a unique  scaling  function  solution  of  (6),  whose  support 
is  [0,  m + 1],  such  that 


2>(*-t)  = l,  x £ 1R.  (9) 

i€7Z 

Moreover,  the  functions  {<p(-  —i),i  € K}  are  linearly  independent  and  totally 
positive  on  IR. 

The  aim  of  this  paper  is  to  construct  new  classes  of  blending  B-bases, 
from  a given  system  {p(-  —i),i  € Zj,  where  Z is  a finite  subset  of  7Z  and  ip 
is  a scaling  function.  In  Section  2 some  preliminaries  are  outlined,  whereas  in 
Section  3 this  construction  is  specialized  to  the  new  classes  of  scaling  functions 
introduced  in  [10].  Finally,  Section  4 is  devoted  to  some  examples. 

§2.  Preliminaries 

Let  I = [a,/?],  with  a,  fJ  integers,  be  a finite  interval  of  IR  and  let  ip  be  a 
compactly  supported  scaling  function,  whose  support  is  [0,  L\,  associated  with 
a mask  a enjoying  the  properties  i)  and  ii)  of  the  previous  section.  Then,  the 
system  of  n = /?  — a + L — 2 functions 

$ :=  {p(x  — i),a  — L + 1 < i < (3  — 1},  x € [a,/?],  (10) 

constitutes  a blending  (cf.  (9))  TP  basis  in  the  space  {/$  generated  by  itself, 
and  fulfils  some  interesting  shape  preserving  properties. 

Indeed,  because  of  the  properties  of  ip  mentioned  above,  the  basis  $ 
satisfies  the  CH  and  the  VD  properties.  Thus,  $ preserves  monotonicity  and 
convexity,  that  is,  any  straight  line  cuts  the  curve  7$  no  more  often  than 
it  cuts  the  control  polygon  [7].  Further  shape  preserving  properties  can  be 
deduced  by  the  generalized  VD  property  for  TP  bases  (see  [2]). 

It  is  rather  natural  to  wonder  whether  $ is  a B-basis,  too.  To  this  end 
we  can  use  the  following  proposition  from  [4], 
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Proposition  A.  A TP  basis  B — (Co,  ...,£„)  is  a B-basis  if  and  only  if  the 
following  conditions  hold: 


inf 


Ci(z) 

0(x) 


x € /,  Cj 


= 0, 


for  all  i / j . 

Clearly,  Proposition  A provides  a useful  test  to  check  if  a TP  basis  is  a 
B-basis.  If  the  check  fails,  one  can  construct  the  unique  blending  B-basis  of 
the  space  [/<j>  by  means  of  the  procedure  given  in  [4,  Th  3.6  and  Th.  4.2]. 


§3.  Construction  of  B-bases  of  Scaling  Functions 

One  of  the  main  advantages  of  the  cardinal  B-spline  as  scaling  function  is 
that  its  mask  has  an  explicit  expression  (cf.  (5)).  A wide  generalization  of 
the  cardinal  B-splines  was  developed  in  [10],  where  a new  family  of  scaling 
functions  has  been  introduced  by  means  of  a new  family  of  masks,  which 
have  an  explicit  expression.  These  scaling  functions  depend  on  certain  free 
parameters,  have  prescribed  smoothness  and,  as  for  the  cardinal  B-splines, 
are  compactly  supported,  totally  positive  and  centrally  symmetric.  They  were 
introduced  as  follows. 

Let  H denote  the  set  of  all  compactly  supported  and  centrally  symmetric 
masks  whose  symbol  is  a Hurwitz  polynomial.  In  [10]  it  was  proved  that  a 
mask  a belongs  to  H if  and  only  if  its  coefficients  are  of  the  type 


(m.fc) 


i = 0, 1, . ..  ,m  + 1,  (11) 


where  m = 2, 3, . . .,  k is  an  even  integer  such  that  1 < k < m,  and 


ik(r)  _ t(r-i ) _ A 2r  + 2\  (r_!) 

* 4 \i-r  + 1 Jr-1  ’ 


r = 0, 1, . . . , A,  K:=k- 
i = r + 1,.. . , K + 1, 


1, 


(12) 


and  b\°\ 


i = 0, . . . , k,  are  such  that 


b(,0)  =b<°) 

uk-r  ur  i 

r = 0,1,...,  A:, 

K 

6(0)  = 2 k-m 

i=0 

det  (4“-jib. 

i = i, ■•■.?)  > o, 

p=l,...,k 


(13) 


(assume  (])  = 0 for  i < 0 or  i > l). 
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Due  to  the  properties  of  a € if,  the  scaling  function  <pmik,  which  is  the 
solution  of  the  scaling  equation 

771+1 

mW  = aiW,fcVm, k(2x  — *)»  x e (14) 

t=0 

is  compactly  supported  on  [0,  m+1]  and  centrally  symmetric,  and  the  functions 
{Tm,k{-  — i), i € ZZ}  are  linearly  independent,  normalized  and  TP.  Moreover, 
recalling  that  a scaling  function  belongs  to  Cr(IR)  if  and  only  if  the  symbol 
can  be  factored  as 


p(z)  = (z  + 1 )r+1qm-r{z),  qm-r(  1)  = 2 r,  (15) 

(see  [8]),  one  can  prove  that  ipmtk  £ Cm_,fc(IR). 

Remark.  Choosing  suitably  the  coefficients  b[°\  the  <pm,k  reduces  to  the 
cardinal  B-spline  of  degree  m,  and  the  <pm,k  can  be  viewed  as  a generalization 
of  the  cardinal  B-splines.  In  particular,  for  k = 1,  the  unique  family  of  scaling 
functions  that  we  obtain  are  the  cardinal  B-splines.  Moreover,  in  the  case 
when  m = 3,  the  coefficients  of  the  mask  (11)  are  a subset  of  those  of  the 
filters  exploited  by  Burt  and  Adelson  in  vision  analysis  [1] . 

Following  the  procedure  outlined  in  the  previous  section,  any  of  the  scal- 
ing functions  ym,k  can  be  used  to  construct  a blending  TP  basis  $m,k  defined 
on  a finite  interval.  Observe  that  a space  is  suitable  for  design  purposes  if  it 
has  a blending  TP  basis. 

By  means  of  Proposition  A,  it  is  easy  to  show  that  the  basis  $m,k  is 
not  a B-basis.  Then  to  obtain  a blending  B-basis  starting  from  the  functions 
^pm,k(x  — i),  we  have  to  apply  the  procedure  given  in  [4],  The  corresponding 
algorithm  can  be  illustrated  as  follows.  Let 

uf  = Pm,k(x  - i),  * = « — to,  — 1, 


where  the  values  of  <pm,k  can  be  evaluated  by  means  of  the  cascade  algorithm 
[12].  For  j = 0, . . . , m — 2,  define  iteratively 


«i+1  :=  < 


u :■  — inf 


H—l 


uj-v  i = m,m  — 1, ...  ,j  + 1, 


Then,  let 


= u?~\ 


i = a — m, . . . ,/?  — 1, 
and  for  j = 0, . . . , m — 2 define  iteratively 


J +1 

. — 


vi  - inf  ( 

^ vi  \ j 

j ) vi+ 1> 

KvUiJ 

* = o,i,.. 

•,/3-2-j, 

.4, 

i = P-  1- 

T- 1 
1 

i 
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The  system  4*®^  = {6;  :=  w™-1,*  = a — m, — 1},  forms  a B-basis.  The 
system  {dibi,i  = a — m, . . . , p — 1},  where  di,  i = a — m, ...,/?  — 1 are  positive 

constants  such  that  da-mba-m  H f dp-\bp-\  = 1,  is  the  required  blending 

B-basis. 

We  remark  that  one  of  the  difficulties  in  applying  this  method  lies  in  the 
evaluation  of  inf(u^/ti;_1)  and  inf(u|/v;+1).  For  instance,  in  the  examples  of 
Section  4,  the  infimums  has  been  evaluated  by  extrapolating  the  values  that 
the  involved  functions  Ui  and  vt  assume  in  a suitable  right  neighbourhood  of 
a and  in  a suitable  left  neighbourhood  of  P,  respectively. 


§4.  Examples 

For  k = 2,  the  mask  (11)  depends  on  a free  parameter  bg°\  which  for  com- 
putational convenience  we  chose  as  a dyadic  fraction:  b^  = 2~h.  Thus,  the 
explicit  expression  of  the  mask  coefficients  becomes 


(j  = 0, 1, . . . , m + 1,  m > 2,  h > m,  — 1),  which  corresponds  to  the  symbol 

= 2-ft(l  + z)m~\z2  + (2h~m+2  - 2 )z  + 1).  (17) 

Observe  that  the  second  term  in  the  mask  (16)  can  be  seen  as  a perturbation 
of  the  mask  of  the  cardinal  B-spline  to  which  (16)  reduces  when  h = m. 

Given  the  interval  I = [a,0],  we  can  construct  the  family  of  blending  TP 
bases 

$m,h  = {<Pm,h(x  ~ *)»«  “ 171  < * < P - !}.  (18) 

where  m > 2 and  h > m — 1.  In  Fig.  1 the  basis  $3  4 defined  on  the  interval 
[0,4]  is  displayed  (dashed  line)  together  with  the  corresponding  blending  B- 
basis  (solid  line)  obtained  by  means  of  the  procedure  outlined  in  the  previous 
section. 

For  k = 4,  the  symbol  p(z ) depends  on  two  free  parameters,  that  is  b^ 
and  b[°\  which  again,  for  computational  convenience,  we  choose  as  dyadic 
fractions:  bg  = 2~h,  b[":>  = 2l~h;  h,l  6 1R  are  arbitrary  numbers  such  that 

h > m — 2 + log2(l  + 2,_1),  in  order  to  fulfil  the  third  of  (13).  Thus,  the 
symbol  has  the  form 

Pm,hAz)  = 2^(1  + z)m~3  (z4  + 2 lz3  + (2~m+4+h  - 2 - 2 ,+1)z2  + 2lz  + l) , 

(19) 

where  m > 3,  and  the  coefficients  aj^\  0 < i < m + 1,  of  the  corresponding 
mask  are 


Also  in  this  case,  the  mask  of  the  cardinal  B-spline  Nm  can  be  obtained 
setting  suitably  the  parameters  h and  f,  that  is,  h = m and  1 = 2.  In  Fig.  2 


the  blending  TP  basis  $5,6,2  defined  on  the  interval  [0, 6]  is  displayed  (dashed 
line)  together  with  the  corresponding  blending  B-basis  (solid  line)  obtained. 
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Fig.  1.  The  blending  B-basis  (solid  line)  and  the  blending  TP  basis  $3,4 
(dashed  line)  in  the  interval  [0,4], 


Fig.  2.  The  blending  B-basis  $^6,2  (solid  line)  and  blending  TP  basis  $5,6,2 
(dashed  line)  in  the  interval  [0,6]. 


Remark.  When  the  scaling  function  is  just  Nm,  the  procedure  outlined  here 
gives  the  basis  of  the  cardinal  B-splines  as  defined  in  [11], 
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